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^ Abstract 

O We further consider a probe fermion in a dyonic black hole background in anti-de 

Sitter spacetime, at zero temperature, comparing and contrasting two distinct classes 
. . of solution that have previously appeared in the literature. Each class has members 

labeled by an integer n, corresponding to the nth Landau level for the fermion. Our 
^ interest is the study of the spectral function of the fermion, interpreting poles in it 

as indicative of quasiparticles associated with the edge of a Fermi surface in the holo- 
graphically dual strongly coupled theory in a background magnetic field H at finite 
chemical potential. Using both analytical and numerical methods, we explicitly show 
how one class of solutions naturally leads to an infinite family of quasiparticle peaks, 
signaling the presence of a Fermi surface for each level n. We present some of the prop- 
erties of these peaks, which fall into a well behaved pattern at large n, extracting the 
scaling of Fermi energy with n and H, as well as the dispersion of the quasiparticles. 
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1 Introduction 



Recent work has considerably enriched our detailed knowledge of the kinds of strongly cou- 
pled physics that can be studied using holographic studies of charged black holes. Such 
studies are in fact quite old, since it was realized very shortly after the formulation of the 
AdS/CFT correspondence [H [3] that electrically charged black holes allow for the study 
of non-trivial strongly coupled physics at finite density or chemical potentialjU E] at both 
non-zero and zero temperature. A rich phase structure was uncovered in refs.|ll E]; fo^' 
(finite volume) systems dual to asymptotically AdS spacetime backgrounds with Reissner- 
Nordstrom black holes in diverse dimensions. (See also the related work in refs.[Bl [7] on the 
phase structure of certain scalar charged black holes.) 

A decade after those studies a new wave of interest emerged, principally due to key ob- 
servations about how charged black hole solutions to asymptotically AdS Einstein-Maxwell 
systems (sometimes coupled to scalars and fermions) may pertain directly to holographic 
models of superconductivity [H IH [10] and to holographic access to Fermi surfaces [TT| [12] . 
(Note that the recent focus has been infinite volume system, therefore focusing on flat black 
holes.) 

In the latter case (refs. pT] [T2] ) the system that displays Fermi surface behaviour 
is a zero temperature Reissner-Nordstrom black hole with a probe fermion. The spectral 
function Gr{u, k) of the probe fermion in this background was shown, at T = and finite 
chemical potential fi, to exhibit a sharp peak at a; = and finite k{= kp). This peak has 
been interpreted as the quasiparticle pole at the edge of a Fermi surface. Here uj is identified 
with Fermi energy above the chemical potential, uj = Ep — fJ^- The dispersion associated 
with this peak was observed to be distinctly non-Landau in a manner consistent with the 
quasiparticles being not effectively free (as in standard Landau theory), but in a strongly 
coupled phase. This is encouraging, since this type of behaviour is reminiscent of that 
which one might hope to gain understanding of in order to gain insight into various novel 
phenomena known from experiment in condensed matter physics (see for example ref . [13] 
for a review). 

Several concerns arise here. Not the least is that the while the results are suggestive, 
a quasiparticle peak alone (especially for a minimally coupled non-backreacting fermion) is 
not damning evidence for a Fermi surface. Some encouragement is to be found in ref. 
where they showed that at T 7^ the peak broadened as a function of T in the expected 
manner, but much more evidence is desirable. Another concern is that the parallel studies 
of holographic superconductivity in the literature all strongly indicate that a charged scalar 
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profile (representing the condensate/order parameter) would have developed at a critical 
temperature well above T = 0, and so an extremal Reissner-Nordstrom solution is presum- 
ably the wrong object to be studying for the T = physics at /i 7^ 0. (In fact, that an 
extremal black hole is probably excluded from this point of the phase diagram was already 
suggested a decade ago[3], since there is an expected decay channel due to super-radiance 
effects [15], inherited from the fact that the charged black holes are really dimensional reduc- 
tions of spinning D3-branes.) The approximate holographic superconductivity studied in 
refs. [HI El [ID] suggest that at T = the system would be (at best) a scalar-profiled charged 
black hole {i.e. not Reissner-Nordstrom) and so one might argue that this is the system to 
which the spectral analysis should be applied. However, further study in a fully backreacted 
M-theory context [161 113 IlHl [19] has shown that the black hole nature of the system entirely 
vanishes at T = 0. An AdS background re-emerges at T = instead. This may all fit well 
with the expectation that if the system has become superconducting below some T = Tc 
the Fermi surface may well no longer be available, although it remains to be seen what the 
appropriate study might reveal. 

With all these concerns up front, it is our view that the work of refs. [HI [12] remains 
of considerable value. It is a simple, remarkable system with many of the key features of 
interest for some of the studies, since away from T = 0, at high enough temperature, a 
charged black hole will still be relevant. The techniques developed and insights gained in 
this study, even if the system is not exactly the right one, will be of value in the broader 
context. 

In this spirit, we continue the study that we began of the Reissner-Nordstrom system 
and the putative Fermi surface in the presence of a background magnetic field "H. In our 
first paper [20], we studied a dyonic black hole, the magnetic component now playing the role 
of an external magnetic field of the system. On general grounds, one expects the system 
to develop an infinite set of Landau levels, the quantized energy levels of the fermions in 
the magnetic field. Our first paper studied the lowest Landau level, and already we saw 
interesting physics, including multiple quasiparticle peaks and a magnetic field dependence 
in the dispersion characteristic of the peaks. 

The non-zero magnetic field case has peaks appearing at w = Ep — fi ^ since 
now the energy of the fermions has been increased as compared to the "H = system, for 
a given /i. So even for the lowest Landau level, the peak is at a; 7^ 0, as we saw in ref.[20]. 
This lowest Landau level has a boundary condition for the probe fermion at the horizon that 
is constructed from the zeroth Hermite function. In this paper we will study more general 
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boundary conditions corresponding to the nth Hermite function, and we find a family of 
quasiparticle peaks signaUing the presence of a Fermi surface at each such n. Each value 
of n, for fixed chemical potential /i, is a distinct Landau level for the probe fermion. 

There are two distinct classes of solution that we present and study here. One class is 



separable (section 3.1) and the other can be understood as an infinite sum of the separable 
solutions (section 3.2). The first, adopted in the work of refs.[2T|l22]. does not have a smooth 
limit to the H = case (for general momentum k^). As a result, we find it less compelling 
as a basis for the generalization of the results and methods of ref. [TT] for the purposes of 
identifying Fermi surfaces by searching for quasiparticle peaks, since the method requires that 
the location (in (w, k)) of the peaks in the spectral function Gr{uj, k) (generically appearing 
at non-zero kx) be read off from the system, and not fixed by hand. We stress that there 
is nothing wrong with studying the separable class of solutions, for the appropriate kind of 
physics of interest. Rather, we are noting that it does not connect to the = physics in a 
way that allows an analysis of the Fermi surface physics in the spirit of the prototype case of 
ref.[TT]. Following from this, we note in section 3.1 that the choices made in ref.[5T], when 
taken to their logical conclusion (performing the required Fourier transform to construct the 
spectral function), lead to no quasiparticle peaks at all. By contrast, the infinite-sum class, 
which we began the study of in ref. |20J, has a smooth limit to the H = case, as we discuss 
in section 3.2, and we carry out the search for quasiparticle peaks and display our results in 
section HI 



2 Free Fermions in a Magnetic Field 

To set the stage and for later comparison with the more complicated case, we begin by 
studying a free fermion with mass m charged under a global U{1). The Dirac action is: 

SD = i j d^x {^YV^^^ - m^^) , (1) 

where "D^ = — iqA^ is the covariant derivative. The classical equations of motion are 
given by: 

i-ff'V^ - m) ^ = . (2) 

Let us first consider the case of zero magnetic field, taking only At = —fi to be non-zero. If 
we choose: 

7° = -ias , 7^ = 0-2 , 7^ = -ai , (3) 
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and take an ansatz of the form: 

then the equations of motion reduce to (for m = 0): 

{uj - qji) (j)± - k^(j)^ = . (5) 

The solutions are then given by: 

kx — ou — q/i with 0+ = 0_ , (6) 
—kx — UJ — q/i with 0+ = — 0_ . (7) 

Note in particular that ai k^ = 0, the non-trivial solution is simply uj — qfi = Q with no 
restriction on 0±(w, k^). 

Let us now consider turning on a magnetic field. Taking A^. = —yH, and following 
a similar procedure as above (with 4>± now depending on y as well), we find as equations of 
motion: 

dy(j)± + m0=p ± {kx + Hqy) (p^ = {uj - qfj,) 0^ . (8) 
Let us assume at present that Hq > and m — and define a new coordinate r] as: 

(9) 



V=VHq(y+^^ , 



in terms of which the equation of motion simplifies to: 

(a,0± ± rjcP^) = ± (a; - q^i) 0^ . (10) 
If we consider solutions of the form: 

0+ = c{uj, kx)In{v) : (t>- = c(a;, kx)In-i{v) : (11) 

where /„ are the standard Hermite functions defined in terms of the Hermite polynomials 
Hn as : 

IM = N^H^ = ^L^ e-^"/^HM , (12) 
we have the requirement that: 

\/2Hqn — uj — q/i . (13) 

Another solution is given by: 

(f)+ = c{uj,kx)In{ri) , (f>- = -c{uj,kx)I„-i{ri) , (14) 



The solutions for the fields are the same as before, but now we have the requirement: 

— ■\/2Hqn = u — qfi . (15) 

Let us now consider the case of Hq < 0. The analysis is mostly the same, except that we 
now define rj as: 



^ = -y/-Hq[y + j^] (16) 



This modifies the equation of motion to: 



V-Hq idr,(p± =F #±) = =F (w - g/i) 0=f • (17) 
The solutions to this equation are given by: 

(f)+{u,k^,y) = c{u,k^)In-i{r]) , 4>-{^,h,y) = c{u, k^)In{r]) , 

(18) 

^—2Hqn = u — ^q , 
and 

0+(cu, k^, y) = c{u, k^)In-i{r]) , 0_(a;, k^, y) = -c{u, k^)In{rj) , 

(19) 

—y^—2Hqn = u — fiq . 

In particular, we note that fiipping the sign on the magnetic field interchanges the r] depen- 
dence for 0-1- . 

Finally, let us study the zero magnetic field limit, which is a delicate issue. We first 
point out that the coordinate change presented in equation ^ is singular in the limit of zero 
magnetic field, so it is not a particularly good basis for studying this limit. However, one 
might think that in order for the restriction on u to match, we have to choose: 

hm^ (±y/2\Hq\n^ = ±k^ = fixed . (20) 

However, in this limit, the functions /„ behave as: 

Iniv) ~ 7r-h"-3^'3e^*=^^-? , (21) 

so both (j)^ and 0_ vanish in the limit of large n. This does not match the zero magnetic 
field solution we found earlier. We can perhaps compensate for this vanishing by including 
an n-dependent overall normalization of our fields but we emphasize that this issue is 
not present when k^ = 0. There, no special limit needs to be taken on n and the Hermite 
functions go to a constant in the limit of — )■ 0. This sort of property will occur in our 
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full probe fermion problem later, for one class of solutions (separable) that we will discuss. 
There will be another class of solutions (infinite-sum) that will not suffer from this limitation, 
having a smooth H = limit for arbitrary k^. This latter class is naturally amenable to 
being used to search for Fermi surfaces by looking for quasiparticle peaks in the spectral 
function, since these peaks will naturally occur at non-zero which we wish to read off as 
output, and not fix a priori. 



3 Probing the Black Hole 

Our probe is a Dirac fermion in an asymptotically AdS4 dyonic black hole background. It is 
charged under the background's U{1), and the metric and fields are given by: 

ds' = '^{-f{z)de + dx' + dy') + '^^ (22) 

F = 2Ha^dx A dy + 2Qadz A dt , 
f{z) = l+{H' + Q')z'-{l + H' + Q')z' = {l-z){z' + z + l-{H' + Q')z') . 

The parameter a has dimensions of inverse length. The negative cosmological constant sets 
the length scale L, A = — S/L^, and the Einstein-Maxwell action is: 



(23) 



where we use signature ( — h ++), and kI = SttGn- The mass per unit volume and temper- 
ature of the background are given by: 

e = —[l + Q' + H'], T = -[3-{Q' + H')]. (24) 

Hi; 4:71 



The coordinates we have chosen above in equation (22) are such that the boundary of AdS4 
(where the dual (2 + 1) -dimensional theory is defined in the ultraviolet) is at z = while 
the horizon of the black hole is at z = 1. 
We choose a gauge such that: 

At = 2Qa{z - 1) , and = -2Ha^y , (25) 

which sets a chemical potential fi = —Qa and a magnetic field "H = —2Ha^. We choose to 
work at zero temperature, which restricts Q and H to satisfy: 

H^ + Q^ = 3 . (26) 
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As discussed in our first paper [20], it is important to realize that this equation does not place 
a restriction on the magnetic field for a given value of Q. The important quantity to hold 
fixed while varying the magnetic field is the chemical potential /i, and as Q decreases for 
increasing H, a increases to compensate, resulting in the physical magnetic field "H being 
able to run its full natural range from zero to infinity. 

Now z is dimensionless in our equations above while all the other coordinates are 
dimensionful. For convenience, in what follows we will rescale our fields and coordinates to 
dimensionless quantities: 

t ^ t/a , X — )■ x/a , y ^ y/a , At ^ aAt , and — )■ a'^A^ . (27) 

The Dirac action for our probe fermion is: 

SD = ij d'^xV^ (^T^^Dm^ - m^^) , (28) 

with the covariant derivative Dm given by: 

Vm = Om + ^cOabAiT''' - iqAM , (29) 

where we have used vielbeins e^^ to exchange curved spacetime indices {M, A^} for tangent 
space indices {a, 6}, and 

= i[r", r''] , Ua^M = e^dMem - eaJve^F^^ . (30) 

The fermion couples to an operator of dimension A = m + d/2 in the (2 + l)-dimensional 
theory. We will later choose the case m = for much of the treatment we will do in this 
paper. We choose for our Gamma matrices: 
















-ia3 \ 




J 



-(T3 
ag 

where the cxj are the standard Pauli matrices. 

We write the Dirac spinor in terms of two 2-component spinors: 

and the equation of motion reduces to: 

^ (d^(f)i + + ^uascpi ± fxa (dy(f)2 + {k^ + 2Hqy) 0i j = , (33) 



with: 



u 



^(cu + 2qQ{z-l)) , 

-gtt V 9zz 



-9tt y/7 

Near the AdS boundary, the solutions to this equation asymptote to: 



\im(j)i{y,z) = ai(y) 

2->0 ^ 2 



We can define 4-component spinors (f)± that have eigenvalues ±1 under 



(34) 



(35) 



(36) 



In the limit of z — > 0, (j)± asymptote to: 

/ -bi \ 
hi 



(l>+iy, z) = z'^ 



-62 

V &2/ 



0-(2/, z) = z" 



ai 

V «2 / 



(37) 



Using the prescription of ref. |25 to calculate the retarded Green function (the spectral 
function of interest), we write: 
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where a and h are the Fourier transform in the y-direction of a and h respectively and hence 
find for the retarded Green function: 
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(39) 



The dual (boundary) field theory should have half the components that the bulk theory has, 
and this is reflected in the Green function by having only two independent terms. Finally, if 
we write: 



5i 



and 



^2 
B2 



then the functions (ai , h\ ) can be written: 

^ 2' -i' 

% (y) = Jim Y (y, e) + Ai {y, e) j , 61 (y) = lim — {y, e) - Ai (y. 



(40) 



(41) 



(with analogous expressions for the Fourier transforms (a, b), with {A, B) instead of (A, B)) 
such that the relevant Green function quantities are simply: 

G^^\uj, K, ky) = limze-''" ^'"^' , Gg^(w, ky) = limze'''" ^' ~ . (42) 
•=^0 Bi + Ai "^^o B2 + A2 

To solve for the fields Ai.Bi, we define: 

2' 2' 

The equations of motion for the fields {A±, B±) when m = is given by: 



Qxx 
9zz 



(^^ ^ Q{l-zy ^ ' 18^(?- zT ) " -^MA+ + ^(9,5_ + (2i^gy + A;,.)S+) 

+ ^ + - -..A.-.(.,i., + (2^,. + Mi^_) 

5+ = - i {dyA_ + (2i7gy + k^)A+) 

B_ = +iuB_+i{dyA+ + {2Hqy + k^)A_) 



Qxx 
9zz 



Qxx 
Qzz 





iuj 


+ i- 


-QqQ + 4uj 


6(1 -z)2 


18(1 -z) 




iu 


+ i- 


-QqQ + 4uj 


6(1 -z)2 


18(1 -z) 




iu 


+ i- 


-QqQ + Aio 


6(1-^)2 




18(1 - z) 




iu 


+ i- 


-QqQ + Aio 


6(1 -z)2 




18(1-2) 



(44) 



These are the equations of motion we presented in ref. ^Dj- Expanding the equations of 
motion at the event horizon, we find (for 7^ 0) the following conditions: 

A±{y,l) = 0, (45) 
\/Q 

d,A^iy,l) = T^ik,B^iy,l) + dyB^iy,l) + 2HqyB^iy,l)) , (46) 

d,B±{y,l) = ^{Tl8VQ{k,d,A±{y,l) + d,{dyA^{y,l) + 2HqyA±{y,l))) 

+ {A8qQ-23u)B±{y,l)) . (47) 

(There are analogous equations for u = 0, but we do not list them here since we will not find 
quasiparticle peaks at a; = 0.) There is some freedom in choosing the boundary conditions 
for the fields, leading to very different physics. We present two particular choices in the next 
two sections. 

3.1 Separable Solutions 

We restrict ourselves at present to qH > 0. Changing coordinates to: 



(4J 
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we can write the y-dependent parts of our equations above as: 

dyB± + {2Hqy + h)B^ = ,/2Hq (d.B^ + r]B^) , (49) 
dyA^ + i2Hqy + K)A^ = ,/2Hq idr,A^ + r]A^) . (50) 

There are two possible ansatze that we can make: 

ansatz 1 : A^ = -A+ = -YA{y)ZA{z) , B. = B+ = Yb{v)Zb{z) , (51) 

or 

ansatz 2 : A. = A+ = YA{y)ZA{z) , 5_ = -B+ = -YB{y)ZB{z) . (52) 

We emphasize that the ansatz with A_ = A^ and B_ = B^ is not consistent with the first 
order equations. Let us start with the first ansatz. If we write: 



dr,YB + vYb = V2nYA , (53) 
dr^YA-vYA = -V2^Yb, (54) 

then the solution takes the form: 

Ya = In-M , Yb = Iniv) , (55) 

where /„ is the Hermite function defined in section [2] Note that /_i = 0. With this ansat^ 
the equations of motion in the bulk become: 

^^ + ^71 ^2+^^^71 v)^^ = -luZA + 2l^JHqnZB , 

Qzz V 6(1 - zy 18(1 - z) J 

~ \ + ^//^ ^o + '^-T^7T~~V~] = +iuZB - 2i^jHqnZA ■ (56) 
9zz \ 6(1 - zY 18(1 - z) J 

If we consider the second ansatz, we can write: 

d,,YB-riYB = -v^Ia, (57) 
driYA + vYa = V2^Yb , (58) 



which has solutions given by: 

YA{y) = In{v), YB{y)=In-M, (59) 



^These are the separable solutions discussed in refs. [211 122], and they are somewhat analogous to the 
standard separable free fermion case reviewed in section [2] We considered them in our work in ref . [20^ for 
n — 0, but focused on the infinite-sum solutions for the rest of that paper. The lack of a smooth H — 
limit for =/= 0, combined with the desire to seek quasiparticle peaks at non-zero k^, led us away from these 
separable cases. We will discuss this more below. 
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and the equations of motion reduce to equations ( 56 ) , the same set of equations of motion that 
resulted from the first ansatz. So in this case, the behavior in the z direction is independent 
of our choice. We will return to the issue of this freedom later. 

We proceed with calculating the Green function given by equation (42), where all the 
fields have been Fourier transformed in both [x^y) coordinates. Note that in the analysis of 
ref. [21], no Fourier transform was performed for the y coordinate. The Fourier transform of 
the Hermite functions is given by: 

dye-^^^yiM = e^^^N^{-zTe<^)^''H^ {-^^ , (60) 



where the normalization Nn was given in equation (12). 

Let us compare to the results and methods of ref. [TT], where H = 0. There, the 
choice ky = is made, and we can examine the physics of the solution we have here at 
H ^ for ky = to see how it connects. If we restrict to ky = 0, we find that for n odd, the 
Fourier transform of In{v) vanishes. Therefore, when n is even, only B±{kx, ky = 0) survive 
when using the first ansatz, whereas only A±{kx,ky = 0) survive when using the second 
ansatz. When instead n is odd, only A±{krc,ky = 0) survive for the first ansatz and only 
B±{kx,ky = 0) survive when using the second ansatz. Therefore, if we wish our retarded 
Green function to be positive, we must use the first ansatz to describe even modes and the 
second ansatz to describe odd modes. This gives: 

GK(c^,A;,,fc, = 0)W'(2)=^. (61) 

Let us now consider the H limit of our solutions. As shown earlier for the free fermion 
in 2+1 dimensions, the Hermite functions in terms of rj do not have nice behavior as H ^ 0. 
Furthermore, the Green function presented above in the if — )• limit cannot match the 
results presented in ref. [11], since their results at if = are a constant only at k^ = 0. 
Therefore, these results suggest that this choice of solutions (used, crucially, in ref. [21] to 
discuss Fermi surfaces) is not connected to the pole found in ref. [TT] at ii = 0, but they 
form a separate set of solutions present at non-zero magnetic field. 

We complete our discussion of this set of solutions by considering the case of qH < 0. 
The coordinate change is now given by: 

(62) 
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Under the change of variables, the j/-dependent part becomes: 

dyB^ + {2Hqy + k^)B^ = ~-^^2Hq {dr,B± ~ 7]B^) , 

dyA^ + {2Hqy + h,)A^ = -^-2Hq {dr,A±-7]A^) . (63) 

If we now consider our two ansatz from before, we find that ansatz 1 results in: 

B+{y,z) = B_{y,z) = ZB{z)In-i{v) , A+{y,z) = -A_{y,z) = Za{z)IM , (64) 
and ansatz 2 gives: 

B+{y,z) = -B_{y,z) = Zb{z)IM , A+{y,z) = A_{y,z) = ZA{z)In-M . (65) 

Both these ansatze give exactly the same equation of motion as before with the appropriate 
modification under the square root: 

-QqQ + Au^ 



9xx 
9zz 

9zz 



too 
6(1-^)2 



18(1 - z) 
-6qQ + Aco 



—iuZa + 2i\/ —HquZs , 



— 2i\/ —HqnZA 



(66) 



6(1 zf ' ' 18(1 - z) 
In turn this suggests that when Hq < 0, the even modes are associated with ansatz 2 and 
the odd modes are associated with ansatz 1. This is the opposite of what occurred when 
Hq > 0. 



Let us summarize our results so far: 



ansatz l{n even) 



qH>0< 



V 



—ZaIu-i 

ZbIji 



ansatz l{n odd) 



ansatz 2{n odd) 



/ ZaIti 

ZbIu-i 
ZaIu 
\ —ZbIu-i I 



qH <0< 



ansatz 2(n even) 



/ ZaIti \ 

Zsln-i 
—ZaIu 
\ ZbIji-i J 

( ZaIh-i \ 

Zsln 

ZaIu-i 
V -ZbIu I 



(67) 

Note the close similarity between these results and the results found for the free fermion 
in 2 + 1 dimensions. Wc find here that for a given ansatz, the single Dirac spinor in the 
(3 + l)-dimensional background contains both solutions that were found for the free fermion. 
Under the flip of the magnetic field for a given ansatz, wc find the same flipping of the r] 
dependence between each 2-component spinors. In addition, the flipping of the magnetic 
fleld interchanges the roles of the two ansatze for describing the modes, which suggests that 
they may describe aligned vs. anti-aligned couplings. These results indicate that the even 
modes and odd modes should be understood as describing two distinct ladders. 
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3.2 Infinite— Sum Solutions 



The separable solutions presented so far have two unfortunate features. Because of the 
solutions' dependence on the coordinate t], the solutions have singular behavior in the zero 
magnetic field limit. To connect these solutions to the zero magnetic field solution, one 



has to take a limit similar to that of equation (20), where the Landau level label n goes to 
infinity: 

lim ^/AHqn = (68) 

i?— s>0,n— s>oo 

Independent of this fact, considering the separable solutions alone leads to a trivial Green 



function, as presented in equation (61). 



We wish to proceed by studying a class of solutions that avoids these two features. In 
order to avoid singularities associated with the coordinate t] , we wish to consider solutions 
that depend on the coordinate y as opposed to rj. Let us take the case of Hq > 0. Motivated 
by the separable solutions discussed earlier, focusing on the behaviour at the horizon = 1), 
we consider B±{y, 1) = B{y, 1) for even n and = — = B{y, 1) for odd n. For n even 
we take: 

dyB{y, 1) + 2HqyB{y, 1) = ,/AHqnBo{l)In-i{V^y) , (69) 

such that we have B{y,l) = BQ{l)In{\/2Hqy). We can write this ansatz in terms of our 
separable solutions from the previous section. 



B^{y,z) = Y,bt^zP{^)In{v) . (70) 

n=0 

and similarly for the other fields. The coefficients fo^'* can be determined from the field 
behavior at the event horizon. For example, for 1) = lQ{y/2Hqy), the coefficients ft^'' 
are given by: 

^in) ^ e-g^ . (71) 

In the limit of — ?■ 0, the infinite sum remains finite without requiring us to take a limit of 



the form of equation (68). We do not need to take n oo since the infinite sum effectively 
does this by defining a new Landau level labeling in terms of an infinite sum of the Landau 
labels of the separable solutions {i.e. the labeling of In{y/2Hqy) instead of Iniji))- Because 
of the nature of the infinite sum, it is important to first perform the sum and then take the 
limit of if —> to ensure the correct result. We emphasize that the infinite sum gives us 
the freedom to vary the new Landau level label as well as independently while sending 
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— )■ 0, a freedom we did not have with the separable solutions presented earlier. This is an 
important difference from the separable solutions we considered in the previous section. 



By substituting the series form of the fields from equation ( 70 ) into the equations of 



motion, the partial differential equation becomes an infinite series of coupled ordinary dif- 
ferential equations for ^j^"^, where the coupled equations are identical to those of equations 



(56). We have traded the partial differential equation for an infinite number of coupled ordi- 
nary differential equations. We choose to proceed by solving the partial differential equation 
using an ansatz that is not explicitly separable since this is more tractable numerically. 
With this ansatz, we find that: 

d-A±{y. 1) = (kJniV^Hqy) + ^AHqnIn^r{^2Hqy)) Bo{l) . (72) 



We notice that for = 0, we find that dzA±{y, 1) oc In-i{\/2Hqy), which agrees with the 
separable solution from the previous section. Finally, we have: 

+ {48qQ-23u)In) , (73) 

where we have used that: 

dyin ( V2^y) - 2qHyIn ( y^y) = - ^AqH{n + l)4+i . (74) 

Note again that for kx = 0, we recover that dzB±{y,l) oc In{y/2Hqy) which is the single 
separable solution presented earlier. In addition, note that for n = we recover the equations 
that we had in ref. |20]. For n odd, we take: 



dyBiy, 1) - 2HqyB{y, 1) = mqnBo{l)In{^2Hqy) , (75) 



such that we have B{y, 1) = BQ{l)In-i{\/2Hqy). With this ansatz, we find that: 

dzA^iy, 1) = {kxh_,{^2Hqy) + ^m^IrX^2Hqy)) B^{1) , 

dzB± = T-^^^ ( — (klh-i + kx^jAHqnh + kx^AHq{n - 1)1^,-2 + AHquI^-i 

iUo \ CO \ 



The advantage of this choice over the finitely separable choice from the previous section 
is that in the H limit, the equations of motion and fields reduce to those of ref. 
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with no restrictions placed on n. In addition, since at = 0, the fields reduce to the simple 
separable case discussed earlier, we can immediately write the Green function for this special 
case: 

Gr{u, K = 0,ky = 0)(^)'(2) = ^ , WH , (77) 

which connects nicely with the zero magnetic field A; = result of ref . [TT] . Furthermore, for 
kx 7^ 0, ky = 0, the Green function gives a non-trivial result, as opposed to the separable 



case which gave a trivial result (see equation (61)). 



This infinite-sum class of solutions, which we studied in ref. [20] at n = 0, is a more 
natural generalization of the results of ref. [11] and better adapted to the method of searching 
for new quasiparticle peaks, which generically appear at non-zero k^. 

3.3 Comments about Complex cj 

Before proceeding with our numerical results, we would like to make a few comments about 
the variable u. In the usual AdS/CFT dictionary, it is important to recall that u is generally 
complex: 

u = Lo^: — iV (78) 

The quantity F is related to the inverse lifetime of the quasinormal mode and u^, is related 
to the energy of the mode. The imaginary part of u must be negative in order for the mode 
to be stable. In addition, in order for the pole in the Green function to be associated with 
a quasiparticle, the pole must have non-zero residue. In our analysis, we restrict ourselves 
to F = 0, and therefore our search for poles is such that the lifetime of the mode is always 
infinite. Away from the poles, the numerical results should not be taken too seriously since 
away from the poles, one should consider a complex u and not a purely real u. This is relevant 
because in our numerics the poles occur immediately after a region where the imaginary part 
of the Green function is negative, which is not allowed if the theory is unitary. We believe 
that the proper treatment of the complex u in these regions would resolve this. In particular, 
the analysis in ref. [22] shows that the pole for = bounces off the real axis, and in the 
purely real analysis the imaginary part of the Green function remains positive on both sides 
of the pole. Therefore, our result suggests that the bounce does not occur at the real axis 
but somewhere in the positive imaginary u regime. This suggests that for Re(co') less than 
our pole, the mode is unstable. In addition, we have not calculated the residue of our 
poles as this is impossible without having the full complex u description of the quasinormal 
modes. Therefore, our poles have not been shown conclusively to be quasiparticles. However, 
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in ref. [SU], the pole studied there was the direct perturbation of the quasiparticle pole in 
ref. [11] under a magnetic field (meaning that zero magnetic field limit matched) and the 
perturbed pole occurred at non-zero uj. In that case, since the H = pole was a quasiparticle, 
it followed that the H ^ pole would also be a quasiparticle. Therefore, in the analysis that 
is presented next, the fact that the poles occur at non-zero u is not entirely surprising, and 
they have a chance of being real quasiparticle excitations. 



4 Numerical Results for the Quasiparticle Spectrum 



We employ the numerical methods that we presented in our earlier paperj^O] to seek solutions 
and explicitly construct the spectral function GR{u,k), using the more general boundary 
conditions we have discussed in section 13.21 for the infinite-sum class of solutions. We refer 
the reader there for the details of our methods, and here present our observations. 

First, we notice that for the first four n, (n = 0, . . . , 3) we find two peaks (signalling 
the presence of poles in the spectral function) for a given n. We observed the pair at n = in 
our earlier paper [20], and this structure persists for a few levels before disappearing, leaving 
only single poles for every n beyond n = 3. In each case where there is a pair, the poles 
fall into two classes. One class appears to be a deformation of the prototype pole found in 
ref. [11]. This type is found at smaller u. The other, found at larger u and k^, appears 
to be in the same class as the lone poles present at large n, fitting into a smooth curve of 
progression. We present the first class of poles in figure [Tj and the second class of poles in 
figure [2] In all figures, cu* is used to denote the value of to at which the peak is located. (In 
both figures we colour code even n poles with blue and odd n with red.) Notice that for the 
first class of poles the even n poles are at higher u while for the second class of poles it is 
the other way around. The behavior of the first poles (at n = 0) was studied extensively in 
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Figure 1: The first (finite) class of peaks for a specific value of magnetic field set by = —0.1. 
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ref. [20] as a function of H, and we expect similar behaviour here for n = 1,2, 3. Therefore 
we focus on the second class of poles, which persists for all Landau levels n, in what follows. 
We present some results for fixed H and varying n in figure [2] (Recall that the physical 
magnetic field is "H = —2a^H.) We note that the behavior of the poles in vs. n becomes 
more predictable at higher n. In addition, we note that for Hq > 0, we find that the solutions 
with ansatz 2 (associated with the odd ns) have a larger u^, than the solutions with ansatz 1 
(associated with the even ns). 
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Figure 2: The second class of peaks for a specific value of magnetic field set by if = —0.1. (a) The 
dependence on Fermi momentum, (b) The dependence on Landau level number. 



In figure |3] we plot the behavior of the poles at fixed n for varying H, for a few 
sample values of n, choosing of n such that it falls in the regime where behaves regularly 
with n. Using these (and much more) data, we can attempt to fit the behavior of u^, with 
the magnetic field and n, for large n, giving: 

~ 0.71 i70-65^o-63 _ (79) 

For how well this function fits the data, see figure |3j In particular, we find that the fit 
improves for larger H and larger n. The exact numbers associated with the scaling of H 
and n are not as important as the fact that they deviate from the free relativistic fermion 
behavior which would have both scaling with exponent 1/2. 

Another important characteristic of these poles is the rate at which the peak is ap- 
proached as a function of k^. This is the dispersion of the quasiparticles associated with 
that peak, and it can be read off from our numerical searches quite readily. In our earlier 
paper [20], we found that non-zero H modified the behaviour seen in ref.[Trj. In particular, 
a peak at higher had more linear {i.e. Landau-like) dispersion. We observe that the same 
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is true for all n. Note however that from figure |2] it can be seen that, for the second class of 
peaks (the one that persist for all n), at higher n, the peaks are at successively lower values 
of kx (in either the even or odd n sequence). We show some sample dispersion behavior of 
the poles at fixed H but different n in figure |4j 

5 Conclusions 

We have generalized the work of ref. |20] to include an infinite number of excitations labeled 
by n. Different n correspond to different Landau Levels for the probe fermion in the presence 
of magnetic field. We studied in detail two separate classes of solution, leading to distinct 
physics. One class is separable (it is used in the work of refs.[2T| [22] ) and the other is an 
infinite sum of the separable solutions (studied for = in our first paper [20j). Both classes 
are probably physical, but the first class does not allow a smooth H = limit for arbitrary k^, 
and so seems less suitable for discussion of quasiparticle peaks and Fermi surfaces in the spirit 
of ref. [11], which is our focus (and that of ref. [21]). In fact, we argued that the separable 
solutions have a constant Green function ai ky = when the ^/-dependence is included in 
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Figure 4: Dispersion behavior near the pole. The red dot marks the position of the pole. 



the Green function definition, which differs from the argument presented in ref. pTj. 

The infinite-sum solutions have key properties that make them attractive. They limit 
to the H = solution smoothly for arbitrary n, with non-trivial dependence on kx- This is 
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the class of solutions we use to find quasiparticle peaks for all the Landau Levels n. Note 
that since u depends on k^, these Landau levels do not have the usual degeneracy found for 
the free fermion Landau levels. 

We also noted that levels given by even n and odd n form distinct towers, distinguished 
by a relative energy shift (that diminishes as n increases) that suggests a aligned/ ant i-aligned 
coupling to the magnetic field. The difference between even and odd stems from two different 
choices one can make at the event horizon. 

We also noticed that the family of quasiparticle peaks, at sufficiently large n, has a 
dependence on H and n that is different that of the relativistic massless free fermion. In fact 
the exponent lies between the free relativistic value of 1/2 and the non-relativistic value of 
unity. This is consistent with our quasiparticles possibly having an induced mass. 
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